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Abstract 

We describe several techniques developed by the High Resolution Fly’s Eye ex¬ 
periment for measuring aerosol vertical optical depth, aerosol horizontal attenuation 
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length, and aerosol phase function. The techniques are based on measurements of 
side-scattered light generated by a steerable ultraviolet laser and collected by an 
optical detector designed to measure fluorescence light from cosmic-ray air showers. 
We also present a technique to cross-check the aerosol optical depth measurement 
using air showers observed in stereo. These methods can be used by future air 
fluorescence experiments. 
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1 Overview 


Cosmic-ray air fluorescence experiments measure light produced by high-energy 
air showers. The technique is calorimetric, with the atmosphere as the calorime¬ 
ter. The integrated scintillation light produced by an extensive air shower is 
proportional to the the primary particle energy and is almost independent 
of the primary particle composition. However, once this light is produced, 
the amount that reaches the observatory depends on how this light propa¬ 
gates from the shower through the molecular and aerosol components of the 
atmosphere to the light-collection optics of the detector. The beneht of the 
atmosphere that makes the air flnorescence techniqne possible is balanced by 
the challenge to nnderstand how the light propagates throngh the atmosphere. 


Several aspects of the High Resolntion Fly’s Eye (HiRes) experiment address 
this challenge. To rednce atmospheric effects, the experiment is located in a 
remote desert where the average hnmidity and clond cover are relatively low. 
The detector stations are located on hills approximately 100 m above the 
desert floor. This places the optical detectors above many aerosols, inclnding 
low-lying dust and ground fog. To measnre atmospheric effects, the experiment 
inclndes steerable nv lasers that probe the apertnre of the detector while 
the experiment is rnnning. The light scattered ont of a laser beam prodnces 
tracks in the same detectors that measure tracks of uv light from air showers. 
Changes in atmospheric aerosols change the amonnt of light scattered and the 
amonnt of light that reaches the detector. Modeling the detector, the laser, 
and the atmosphere generates a simnlated set of detector measurements. The 
simulated data are compared and ht to the actnal detector measurements of 
the laser. This procedure yields a set of parameters that describe the aerosols. 
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2 Components of the Atmosphere 


Most of the light observed from air showers is produced in the the troposphere. 
With regard to atmospheric fluorescence measurements, the atmosphere has 
two main components: the molecular atmosphere and aerosols. Both compo¬ 
nents can be described in terms of optical depths, denoted tm and ta- Assum¬ 
ing single scattering, the transmission through a vertical column of atmosphere 
can be expressed as T = exp (—tm — ta)- 


2.1 The Molecular Component 


The molecular atmosphere, comprising mainly nitrogen and oxygen molecules, 
has absorption and scattering properties that a re well understood. Be cause 


the 300 to 400 nm wavelengths of interest (see Kahimoto et ah ( 1 99 (tIi and 


Naga.no et ah (|2nn4B are much larger than the nitrogen and oxygen molecules, 
the electric held can be approximated as constant across a molecule, and 
consequently, light prop agation can be described using Rayleigh scat tering 
theory [se e , for e xample. Nans and Ubachs ( 200011 . Miles et ah ( 2001 1. and 
Bncholtz ( I 995 II ]. The molecular component is relatively stable, with small 
variations due to height-dependent changes in temperature and pressure. Be¬ 
cause a purely molecular atmosphere has the the clearest and best-understood 
viewing conditions, it provides a boundary value solution to light-propagation 
studies. Air density can be obtained by applying the ideal gas law to ra¬ 
diosonde measurements (collected at the Salt Lake, UT, and Reno, NV, air¬ 
ports) of pressure and temperature as a function of height. The dependence 
on height follows an adiabatic model, with variations at the 5% level caused 
by temperature and pressure changes associated with different seasons and 
weather systems. For circularily or unpolarized light, the scattering probabil¬ 
ity as a function of scattering angle, 6, follows 1 -f- cos^{6), with higher-order 
corrections that are negligible for the purposes of this work. 


Applying the methods described in Bodhaine et ah (199^ to Salt Lake ra¬ 
diosonde data, yield values for tm at 355 nm, a wavelength that falls near the 
middle of the air fluorescence spectrum, tm from the 1.5 km elevation of HiRes 
(850 g/cm^) through the entire atmosphere is 0.50. This value drops to 0.37 
through 10 km (223 g/cm^) above ground, and to 0.17 through 3.6 km (557 
g/cm^) above ground, with seasonal fluctuations of 0.005. For comparison, we 
note the average ta between ground and a 3.5 km height was measured at 
355nm to be 0.04 (see section 4). 
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2.2 The Aerosol Component 


The aerosol component is more complicated. Aerosol sizes range from 10“® 
cm (large molecules) to 10“^ cm (particulate matter). The distribution of 
aerosols changes with location, height, and time. Local dust can be lifted from 
the desert floor by wind and transported to elevations typically less than 1 
km. More distant sources can create aerosols that may be carried over large 
areas by upper-level air flows. A precise analytic description of absorption and 
scattering is impossible to obtain because of the variable nature of the aerosol 
size distributions and dielectric constants. Light scattering from aerosols of 
a sp ecihc size and spherica l shap e can be described by Mie scattering the¬ 


ory (jBohren and Huffman (jlhSHTl L In general, the scattering probability as a 


function of angle is characterized by a forward peak, relatively little scattering 
at 90 degrees, and some enhancement for backward scattering. However, be¬ 
cause the precise details of aerosol size, shape, and distribution are not known 
a priori, one must rely on phenomenologically based models for guidance. A 
commonly used simulation package, MODTRA N (Moderate Re solution Trans¬ 
mission), was developed by the US Air Force ( Longtin ( 1988ll L The package 


contains aerosol models for different regions including a “US Standard Desert’' 
parameterization. 


In contrast to the molecular atmosphere, aerosols are known to vary on short 
time scales. Weather fronts can load the atmosphere with wind-driven aerosols. 
Rain and snowstorms can flush haze and dust from the atmosphere, often leav¬ 
ing near-molecular viewing conditions. The aerosol density at ground level, the 
vertical distribution of the aerosols, and the scattering and absorption prop¬ 
erties can all change within one night of observing. Atmospheric monitoring 
at HiRes therefore entails measuring the time dependence of the distribution 
of aerosols in the detector aperture and determining their gross properties to 
the level required to reconstruct air showers. 


In this article, we describe the measurement of three aerosol quantities: verti¬ 
cal aerosol optical depth (r^), horizontal aerosol scattering length at detector 
level (HAL), and aerosol phase function, r^, measured from the elevation of 
the detector to an elevation of 3.5km, is the most important quantity of the 
three for several reasons. First, because the held of view of a HiRes detector 
begins at an elevation angle of 3 degrees, and because the HiRes stations are 
on hills, none of the aperture extends to ground level. At a distance of 19 
km, the lower edge of the aperture is 1 km above detector station height. Sec¬ 
ond, ta describes the total aerosol loading, and most of the aerosols lies below 
the height at which most of the light is observed from air showers, especially 
for showers above 10 EeV, which tend to be further away. Therefore, under 
the assumption of horizontal aerosol uniformity, (also addressed in this arti¬ 
cle), most of the light reaching the detector is sensitive to the total amount 
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of aerosols; the correction for aerosol transmission depends more on eleva¬ 
tion angle than on distance. As viewed along the optical axis of a lowest ring 
HiRes detector, an uncertainty of 0.02 in (r^) corresponds to an uncertainty 
of 11% in transmission. Over the full 30-3 degree range of elevation angles, 
this uncertainty changes from 4% at the very top to 32% at the very bot¬ 
tom. For air-showers, these uncertainties translate to an uncertainty in energy 
reconstruction of order 10 to 20%. 


3 Apparatus 


Measurements of ta and HAL are determined from laser tracks recorded by 
the HiRes detectors. 


3.1 The Steerable Laser System 


The steerable laser system, dubbed HR2SLS Wiencke et ah ( lOOOll . hres a 
pattern of 1000-1300 shots per hour. Vertical shots can be used to extract 
vertical aerosol optical depth. Horizontal attenuation lengths can be extracted 
from near-horizontal shots hred across the HiResl detector station. A pattern 
of elevated shots hred in 15-degree steps in azimuth may be used to study 
atmospheric horizontal uniformity. 


The laser produces a 7 ns wide pulse of light at 355 nm. This wavelength 
closely matches the 357 nm huorescence line. The maximum beam energy 
is 7 mJ. The beam can be steered in any direction above the horizon. The 
direction, energy, and time of each pulse are measured and recorded locally. 
The laser bring times are synchronized using a GPS clock so that laser events 
can be isolated from other HiRes triggers in a physics-blind manner. Different 
energies and polarizations can be selected using biter wheels that provide a 
combination of quarter wave-plates and attenuation biters. 

One percent of the beam is sampled by a photo diode detector to obtain a 
relative measurement of the pulse by pulse energy. Special measurements are 
taken at regular intervals to determine the absolute calibration. For these mea¬ 
surements, an energy-absorbing pyroelectric probe is placed temporarily where 
the beam exits to the sky. The pyroelectric probe is periodically recalibrated 
to 5% by the manufacturer using an NIST traceable system. 


Circularly polarized beams with energies from 0.05 to 3 mJ were used in 
obtaining measurements. Circular polarization was chosen over linear polar¬ 
ization to reduce polarization dependent scattering ebects and simplify simu- 
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lations. The energy range must match the detector dynamic range for different 
beam directions. For example, vertical shots needed to measure ta are viewed 
at a distance of 12.7 km. These shots require energy of 1-3 mJ. On the other 
hand, the horizontal shots needed to measure HAL must pass within a few 
hundred meters of the detector. These shots require a much lower energy, typ¬ 
ically 0.1 mJ, to avoid detector saturation. For each combination of hlter and 
direction setting, the laser hres 15 shots which are later averaged to reduce 
statistical uncertainties. 


3.2 The HiRes Detectors 


The detectors that view the laser tracks are described briefly. HiRes consists of 
two stations, HiResl and HiRes2, which are separated by 12.6 km. The analysis 
described here used detectors (Abu-Zayyad et al. 2000) at the HiResl station. 
These detectors are arranged to view nearly 360 degrees in azimuth from an 
elevation angle of 3-17 degrees. A spherical mirror (3.78 m^ effective area) at 
each detector reflects light from the night sky through a uv pass hlter onto 
a cluster of 256 hexagonal photomultiplier tubes (PMTs). Each PMT views 
about 1x1" of night sky. As a laser light pulse passes upward through the held 
of view of a mirror, light scattered from the beam is collected and focused to a 
spot that crosses downward across the cluster, triggering PMTs along its path. 
The nominal spot is about the size of a PMT. The nominal PMT gain is 10^, 
and the nominal quantum efficiency is 25%. The signal in a PMT from laser 
shots (and air showers) is typically hundreds to thousands of photoelectrons. 

The detector readout system at HiResl uses sample-and-hold electronics. Sig¬ 
nal digitization and data formatting are performed by local electronics at each 
mirror. PMT outputs are AC coupled. The signal-processing channel used for 
this work has a three-pole Bessel hlter with a 375ns time constant and a charge 
integration gate of 5.6 ps. TDCs measure times between PMT triggers by inte¬ 
grating a gated constant-current source. Discriminator thresholds are adjusted 
every 4 seconds to maintain a count rate of 200 Hz. The trigger requirement 
for a signal digitization is a sixfold coincidence between PMTs. Typical trigger 
rates are 5-10 Hz per mirror. System dead-time is less than 2%. 


4 Determining Vertical Aerosol Optical Depth 

Under the assumptions of horizontal uniformity, and single scattering, the 
transmission, Ta, through an optical depth of aerosols, ta, at an elevation 
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Fig. 1. Geometry for measuring aerosol optical depth. 


angle 6 is 


Ta 


exp( 


Ta 

sin(0) 


) 


( 1 ) 


This equation includes both absorption and scattering. For the case that all 
the aerosols are between the starting and ending elevation of interest, and 
the aerosols are uniform horizontally, Ta is independent of the aerosol vertical 
distribution. 

The method to determine ta normalizes the measurement of a laser beam to 
the measurement that would be expected under aerosol-free (i.e. molecular) 
atmospheric conditions. Figure 1 illustrates the geometrical configuration of 
the laser and the detector. 

The amount of light from the laser that reaches the detector can be written 
as 


Nobs — N^lTmiTai{Sa + Sm)Tm2Ta2 


( 2 ) 


where Nqbs is the number of photons observed, is the number of photons 
emitted per laser pulse, Tmi and Tai are the molecular and aerosol trans¬ 
mission components from the laser to the point of scatter, Sa is the frac¬ 
tion of photons scattered toward the detector by aerosols, Sm is the fraction 
scattered by the molecular component, and Tm 2 and Ta 2 are the molecu¬ 
lar and aerosol transmission components from the point of scatter to the 
detector. Normalized to the prediction for a purely molecular atmosphere 
{Nmol = N^^TmiSmTm 2 )-, equation 2 becomes 


Nqbs 

Nmol 


TaiTa2{1 + 1 ^) 


( 3 ) 
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Applying equation 1 yields 


^^(Im 

sin 62 


Nobs ^ , Sa.. 

- - )-ln(l + —)) 

^^MOL JM 


(4) 


The two unknowns, ta and Sa, illustrate the ambiguity between scattering 
and transmission inherent in conventional LIDAR measurements. The ambi¬ 
guity can be eliminated if the scattering occurs in a region where Sa is small 
compared to Sm- Equation 4 then simplihes to 

sin9isin62 . Nqbs . . 

sinOi -f sin92 Nmol 


In this case, ta becomes a simple function of geometry and the number of 
observed photons that are normalized to Nmol- Nmol can be determined 
either by Monte Carlo simulation, or from data collected under molecular 
conditions. The latter method does not require an absolute calibration of either 
the laser or the detector. However, if the reference data were collected during a 
period when aerosols were present, the measured value of ta will be lower than 
the actual value. Normalizing to simulated data requires absolute calibration of 
the laser and detector and simulation of light propagation through a molecular 
atmosphere. This analysis was performed with data normalization and with 
Monte Carlo normalization. The results are shown in hgures 2 and 3. 

The calculation of Nqbs averages over sets of 15 vertical laser shots hred 
each hour at two energy settings. A region of PMTs was selected at the top 
of the track. For the data-normalized analysis, the region was 2 PMTs wide 
by 4 PMTs high, corresponding to an elevation of 3.5 km above the laser. 
For the Monte Carlo-normalized analysis, a 2x2 PMT grid was used, with a 
corresponding elevation of 3.8 km. 

Nmol for the Monte Carlo normalized analysis included effects of multiple¬ 
scattering, estimated to be at the 5% level for a molecular atmosphere. To 
reduce computation time, the Monte Carlo simulation exploited the symmetry 
of a vertical laser beam. Nmol for the data normalized analysis used data 
collected on December 12, 2000, UT hour 6. 

The data sample in this analysis extends from September of 1999 to May 2001. 
Figure 2 shows the time evolution of hourly ta measurements obtained by the 
two methods of analysis, and the difference between the measurements. The 
projection of these data onto the ta axis and the correlation between them are 
shown in hgure 3. The systematic difference of 0.009 between the two analyses 
could be caused by the presence of aerosols during the reference “clear period” 
or by the relative calibration between the laser and the detector. For reference, 
0.009 corresponds to a calibration uncertainty of approximately 5%. 
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Fig. 2. Time evolution of ta- The top panel (A) shows the data-normalized mea¬ 
surements. The middle panel (B) shows the Monte Carlo-normalized measurements, 
and the bottom panel (C) shows the difference. 

The asymmetric shape of the ta distributions in hgure 3 are consistent with a 
molecular limit (at ta=0), and a tail on the high side due to aerosol loading. 
The means of the ta distributions are 0.036 and 0.046. Their full widths at 
half maxima are 0.035. 

We have identihed three sources of systematic error that dominate the mea¬ 
surement uncertainty. A 10% uncertainty in the relative calibration between 
the laser and the detector contributes a systematic error of 0.02 to the Monte 
Carlo-normalized ta average. A 5% time dependent variation in the gains of 
the PMTs used in the analysis contribute 0.01 systematic error to both mea¬ 
surements. Finally, the reference aerosol-free normalization data contributes 
an uncertainty of 0.01 ta to the data-normalized measurement. The number of 
sets of laser shots analyzed reduces the statistical error of the mean to 0.001. 
Adding the systematic errors in quadrature yields 
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Fig. 3. Vertical aerosol optical depth measured at 355 nm. The upper left (A) and 
lower right (D) panels show the measurement using the data-normalized and the 
Monte Carlo-normalized techniques. Panels (B) and (C) show the correlation and 
the difference between the measurements. 


Mean ta = 0.046 ± 0.022^^5 ± O.OOl^tat (Monte Carlo-normalized) 

Mean ta = 0.036 ± 0.014*^8 ± O.OOlstat (data-normalized) 

We note that eqnation 1 yields a correction factor of 1.25 for light traveling 
along the optical axis of a HiResl detector through a ta of 0.04. 

To study the vertical distribution of aerosols, the analysis was repeated using 
sets of PMTs corresponding to different elevations. The mean difference be¬ 
tween Ta for 0-3.5 km and 0-2.5 km was found to be less than 0.0025. See figure 
4. The mean ta for 0-1.5 km was 0.008 smaller than mean ta for 0-3.5 km. 
These numbers indicate that the bulk of the aerosols during this period were 
distributed in the lower 1.5 km of the atmosphere. They are also consistent 
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Fig. 4. Correlation (panel B) between aerosol optical depth as measured from z = 
0 to 2.5 km (panel A) and from z=0 to 3.5 km (panel D). This analysis uses data 
normalization. If the average 0.04 ta was uniformly distributed in height, an average 
difference of approximately 0.01 would have been observed. This “uniform height” 
model is ruled out because the difference (panel C) is much smaller. 

with an average aerosol scale height of approximately 1 km. 

The optical-depth technique can also be applied to inclined laser tracks. To 
test horizontal uniformity, we compared optical-depth measurements from two 
laser directions that were symmetric about the line between the laser and 
the HiResl detector station. The pattern of laser shots includes shots of 15 
degrees elevation angle and azimuth angles of ± 75 degrees with respect to 
HiResl. The portions of the tracks used in this measurement are 3 km above 
the ground,and are separated by roughly 20 km. Although the separation is 
less than 20 km over most of the symmetric light paths (the light starts at 
the laser and ends at the HiResl detector), the symmetric measurements still 
probe different parts of the HiRes aperture and are also made by different 
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Fig. 5. Comparison of vertical aerosol optical depth from 15-degree inclined laser 
shots fired at an azimuth angle of 75 degrees to the north (panel A) and 75 degrees 
to the south (panel D) with respect to HiResl. The regions of the tracks used for 
the measurements are separated by roughly 20 km. The correlation (panel B) and 
difference (panel C) provide a test of aerosol horizontal uniformity. 

mirrors of the detector. The left/right differences in ta fall within ± 0.01 (see 
hgure 5) for more than 90% of the data. A gaussian £t to the difference in 
these Ta measurements yields a sigma of 0.008. In addition, the absolute scale 
of the average ta, as measured with inclined tracks, is consistent with the 
average obtained from the analysis of vertical tracks. 


4.1 Horizontal Attenuation Length and Aerosol Phase Function 

The horizontal attenuation length (HAL) is the aerosol attenuation length at 
ground level. HAL provides a normalization of the aerosol density. The phase 
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Fig. 6. Track made by nearly horizontal laser shot, as recorded by the HiResl 
detector. 

function, or normalized differential scattering cross section, of the aerosols is 
needed to correct for Cherenkov light scattered from an air-shower. HAL and 
phase function are measured simultaneously with a ht to a near horizontal 
shot from the HR2SLS. This shot, passing within 450 m of the HiResl site, 
generates a long track in the detector (see hgure 6) that is observed over a 
large range of scattering angles. Assuming horizontal uniformity, an aerosol 
scattering albedo near unity, and known molecular scattering conditions, the 
longitudinal track prohle can be ht to determine HAL and the aerosol phase 
function (see hgure 7). 

The functional form of the aerosol phase function is 

A oc e-B" + Ce“ (6) 

ail 

where 6 is the scattering angle in radians and B,C,D are free parameters 
in the ht. This functional form reproduces the shape of a Mie theory-derived 
phase function except in the very forward scattering region (< 2°), where only 
a small fraction of the total scattered signal appears. 

Figure 8 shows the measured distribution of ground-level scattering coefficients 
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Fig. 7. Two examples of fits for Horizontal Attenuation Length (HAL) and aerosol 
phase function from data collected under two atmospheric conditions. Panels A 
and C show light intensity recorded from near-horizontal laser shots at the HiResl 
detector as a function of scattering angle. Panels B and D show the corresponding 
aerosol phase function obtained from the fit. For comparison, the molecular phase 
function and the Mie theory-derived phase function of the Longtin desert aerosol 
model are also shown. 


(tto = HAL ^). If the very hazy data (cto > 0.1 km are excluded, the average 
value of tto is 0.04 km“^. 


A mean aerosol scale height can be inferred by dividing the mean ta by the 
mean ckq. This calculation yields 1 km. We note this value is consistent with 
the differences between the ta measurements at different heights as described 
previously. 
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Fig. 8. Measurements of oq from near-horizontal laser shots, ao is the inverse of the 
horizontal aerosol scattering length at the height of the detector. Panels A and B 
show the distribution on different scales. The time dependence is shown in panel C. 


5 Checking the Aerosol Optical Depth using Air-Shower Data 


A sample of air-showers observed by both HiResl and HiRes2 was used to 
cross-check the sca le of the atmospheri c laser measurements. The technique, 
detailed elsewhere ( Abassi et ah ( 20051) 1. is based on the principle that multi¬ 
ple eyes viewing a common track segment should reconstruct the same number 
of emitted photons after corrections for geometrical factors and atmospheric 
attenuation. If two detectors view the same shower through different optical 
depths, an error in the atmospheric model can cause an inconsistency between 
the two detectors in the number of emitted photons from the common track 
segment. Although this technique is a useful cross-check, the detection rate of 
air showers is too low to provide precise measurements over a single night. 
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For this study, a sample of 415 stereo events recorded between August 25, 2000, 
and May 23, 2001, was selected. To ensure well-constrained stereo reconstruc¬ 
tion, the opening angle between the shower-detector planes was required to 
be between 8 and 172 degrees. A cut on the minimum viewing angle between 
the shower axis and the mirror axes of 30 degrees at each detector reduced 
asymmetric contributions from scattered Cherenkov light. For each event the 
light-balance An was calculated from the number of photons reconstructed 
by the HiResl and HiRes2 detectors from the common track segment. 


Atv = In 


( NHiResl 

V NHiRes2 


(7) 


Light-balance is plotted as a function of the linear asymmetry A^, 

Ar = rHiRes2 “ ^HiResl 


( 8 ) 


where rniResX is the distance between a detector and the center of the common 
segment. 

Figure 9 shows this analysis for three aerosol models: 1) near molecular ta 
= 0.001; 2) Ta = 0.040 HAL = 25; 3) “Standard Desert” ta = 0.100 HAL = 
12 km. The Standard Desert parametrization leads to an over correction. On 
average, the far detector reconstrncts a “brighter” common segment than does 
the near detector. Conversely, if one assnmes a complete absence of aerosols 
(i.e. molecular atmosphere), a systematic under correction is observed. The 
0.04 Ta 25 km HAL model, obtained from laser measnrements, falls between 
these extremes, yielding the smallest slope of the three. A ta of 0.055 yields a 
flat slope. The systematic error in this method is approximately 0.01. 


6 Conclusions 


We have developed several methods to measure the atmospheric aerosol prop¬ 
erties of optical depth, extinction length, and phase function. These methods 
use side-scattered light from a 355 nm uv laser that produces tracks in the same 
optical fluorescence detectors that measure tracks from extensive air showers. 
We have measured distributions of aerosol optical depth and horizontal ex¬ 
tinction length at the location of the High Resolution Fly’s Eye experiment. 
We found the distribution of aerosols is horizontally uniform over a distance 
of 20 km to less than 0.01 in ta on more than 90% of the nights observed, and 
that most of these aerosols are in the lower 2 km of the atmosphere. Between 
fall 1999 and spring 2001, the optical depth of the aerosol component is on 
average less than 25% of the optical depth of the molecular component. As a 


16 



0 


z 


TA(model) = 0.00 1 

Slope = 0.0252 
±0.0023 


TA(model) = 0.040 
Slope = 0.0077 
±0.0023 


TA(model) = 0.1 00 
Slope = —0.0230 
±0.0027 


-1 


I_L_ I ■ I _I_I_I 

-15 -10 -50 5 10 15 

A, (km) 



-15 -10 -50 5 10 15 

Ar (km) 



-15 -10 -50 5 10 15 

Ar (km) 


Fig. 9. Light balance as determined by applying three different aerosol models to a 
single set of stereo events: top panel: near molecular, ryi=0.001; middle panel: ta = 
0.040 bottom panel: ta = 0.100; (see text). 


cross check, we note that the average ta, as obtained by laser measurements, 
is consistent an average determined by comparing commonly viewed segments 
of extensive air showers observed by the two HiRes fluorescence detectors. 


7 Acknowledgments 


This work is supported by US NSF grants PHY-9321949, PHY-9322298, PHY- 
0998826, PHY-024528, PHY-0305516, PHY-0307098, and by the DOE grant 
FG03-92ER40732. We gratefully acknowledge the contributions of the tech- 


17 




nical staffs at our home institutions. The cooperation of Colonels E. Fischer 
and G. Hatter, the US Army, and Dugway Proving Ground staff is greatly 
appreciated. 


References 

Abassi et ah, in preparation. 

T.Abu-Zayyad et ah, Nucl. Inst. Meth. A450 (2000) 253. 

B. Bodhaine et ah. Journal of Atmospheric and Oceanic Technology, 16 (1999) 
1854. 

C. F. Bohren and D.R. Huffman, “Absoption and Scattering of Light by Small 
Particles,” New York, Wiley Interscience, 1983. 

A. Bucholtz, Applied Optics Vol. 34 No. 15 (1995) 2765. 

Kakimoto et ah, Nucl. Instr. Meth. A372 (1996) 527. 

D. R. Longtin, “A Wind Dependent Desert Aerosol Model: Radiative Proper¬ 
ties,” Air Force Geophysics Laboratory, AFL-TR-88-0112, 1998. 

R.B. Miles, W.R. Lempert, and J.N. Forkey, Meas. Sci. Technol. 12 (2001) 
R33. 

Nagano et ah, Astroparticle Physics, V22 (2004) 235. 

H. Nans and W. Ubachs, Optics Letters Vol. 25 No. 5 (2000) 347. 

L. Wiencke, et al, Proc. SPIE Vol 3818, (1999) 46. 


18 



